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A concrete strategy is presented for generating strong topological insulators in d + d' dimensions 
which have quantized physics in d dimensions. Here, d counts the physical and d' the virtual di¬ 
mensions. It consists of seeking d-dimensional representations of operator algebras which are usually 
defined ind + d' dimensions where topological elements display strong topological invariants. The 
invariants are shown, however, to be fully determined by the physical dimensions, in the sense that 
their measurement can be done at fixed virtual coordinates. We solve the bulk-boundary correspon¬ 
dence and show that the bormdary invariants are also fully determined by the physical coordinates. 
We analyze the virtual Chern insulator in (1 -t l)-dimensions realized in Ref. [Tl and predict quantized 
forces at the edges. We generate a novel topological system in (3 + l)-dimensions, which is predicted 
to have quantized magneto-electric response. 


INTRODUCTION 

In Refs, m 121, Kraus et al put forward a framework 
in which topological insulators can be generated in ar¬ 
bitrary dimensions but with some of these dimensions 
occurring in a parameter space. For example, a Chern 
insulator defined in one physical dimension and one vir¬ 
tual dimension was proposed and realized in the labora¬ 
tory using coupled one-dimensional wave guides, and 
several interesting topological effects were pointed out. 
The paper states at some point that the entire approach 
could work in any number of dimensions, and indeed, a 
Chern insulator of the second kind in 2 -I- 2 dimensions 
was found in Ref. |[2|. The virtual Chern insulator in 
1 -t- 1 dimensions was since observed in several other 
experiments 110. 

There is certainly a common feature between the 1 -i-1- 
and 2-(-2-dimensional systems mentioned above, namely 
the quasi-crystalline structure, but beyond that they 
seem to be exceptional isolated topological examples of 
this type. One interesting question is if these systems 
actually fit into a larger class unified by a common prin¬ 
ciple? The defining and unifying property was already 
stated in in for (1 -(-1) dimensions, and could be general¬ 
ized as follows: A virtual (d-fd')-dimensional topological 
insulator is characterized by a strong topological invari¬ 
ant in {d + d' )-dimensions but the invariant is measurable 
or computable entirely from the physical d dimensions. 
How can such systems exist? As we shall see, the reason 
is two-fold: 

1. Operator algebras which are naturally defined or 
represented in higher dimensions sometime accept 
faithful {i.e. no information is lost) representations 
in lower dimensions. 

2. A certain self-averaging property occurs in these 
lower dimensions, which has the effect of integrat¬ 
ing out the virtual dimensions. 

These principles were already stated in |[lll2|, but since 
there was no explicit expression of the invariant, there 


was some confusion and, for example, the claims were 
challenged in 0. In their response letter ||7l, the authors 
reiterated the statements but again no expression for the 
invariants was provided. They conceded in this letter 
that the boundary phenomena depend on the virtual 
dimensions, alluding that the topological class cannot 
be determined from boundary phenomena without ex¬ 
ploring the virtual dimensions. In this work, we use an 
operator algebra formalism to provide the missing ex¬ 
pressions of the topological invariants. We pin-point the 
fundamental difference between the irrational and ratio¬ 
nal cases, which is the ergodic vs non-ergodic character 
of the lattice translations. The self-averaging property, 
which exists only for the irrational case, can then be rig¬ 
orously formulated using Birkhoff's ergodic theorem ||8l. 
We also solve the bulk-edge correspondence and prove 
that the boundary invariant, which is equal to the bulk 
invariant, is entirely computable from the physical di¬ 
mensions. In other words, the topological class can be 
determined even if we have access only to the boundary 
of the systems. The formalism also enables us to make 
three new physical predictions for the experimentalists. 

We recall that the bulk-edge correspondence princi¬ 
ple is not just about the emergence of the boundary 
states. The bulk-boundary program was defined and 
solved by Hatsugai 121 for the clean integer quantum 
Hall effect, and consists of 1) defining a boundary in¬ 
variant, 2) proving that the boundary and bulk invari¬ 
ants are equal, 3) determining the physical meaning of 
the boundary invariant. In the general setting, which 
includes disorder and irrational magnetic flux values, 
the bulk-boundary problem was solved by Kellendonk, 
Richter and Schulz-Baldes Eoiini using a K-theoretic 
approach. Our analysis relies on these works, which are 
adapted to the present context without detailed proofs. 
This is mainly because a comprehensive manuscript on 
the bulk-boundary principe for the complex classes is 
in preparation 1121 . For the virtual Chern insulator in 
(1 + l)-dimensions, we find that the boundary invariant 
is equal to the average mechanical force spontaneously 
exerted at the edge of the system. The existence of such 
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quantized forces were first predicted in Refs. IIT3l[T4l . We 
also formulate and analyze the boundary invariant for 
the virtual second-Chem insulator in 3 + 1 dimensions. 

The analysis is carried out at the level of operator al¬ 
gebras, which we find useful for the following reasons: 

• It provides a natural framework for the analysis of 
homogeneous aperiodic systems IITSl . It enables 
us to incorporate incommensurate and random 
potentials yet keep the calculations explicit. As 
such, we are able to translate all the claims made 
in 12 into rigorous mathematical statements and 
to strengthen some of them, especially those re¬ 
ferring to the edge physics. While we agree that 
in condensed matter physics the ultimate confir¬ 
mation comes from experiment, the mathematics 
may be appreciated here because it gives definitive 
answers to some of the open issues. We also feel 
that we can teach the reader in these 13 pages how 
to slightly tweak the old tools and adapt them to 
these more general settings. 

• It provides the big picture, in the sense that seem¬ 
ingly unrelated models can be connected by one 
relatively simple operator algebra. Many operator 
algebras have been studied independently and a 
great deal of information is available (see for ex¬ 
ample IIT^ ). Perhaps the most important piece 
of information is the Ko'group which encodes the 
topology of the projections, in particular, of the 
spectral projections of the Hamiltonians generated 
from the algebra. Then the patterns seen in the 
energy spectrum, such as the Hofstadter butter¬ 
fly, can be qualitatively and quantitatively under¬ 
stood from the Xo'group alone llTl . The fCo'group 
of the rotational algebra, which generates our ex¬ 
ample of (3 + l)-dimensional topological insulator, 
has 8 generators among which one has non-zero 
second-Chem number (equal to 1). As such, even if 
we generate the Hamiltonians randomly, we have 
12.5% chance to generate a second-Chem insula¬ 
tors, hence no numerics are needed to prove that 
such Hamiltonians exist. 

• Lastly, it provides an outstanding computational 
toolbox. Note that no proof (or reference) was pro¬ 
vided in UJ121 for the quantization of the invari¬ 
ants in the strict irrational settings. After all, the 
whole point iTj there was that one cannot "break 
the symmetry" of the quasi-crystal and go into the 
rational setting where one can use a range of ar¬ 
guments. The index theorems for the irrational 
settings, which give the quantization of the first 
IfTSl and higher |[T9l Chem numbers, were estab¬ 
lished using the operator algebra formalism. Also, 
for example, it is well known that one cannot use 
perturbation theory even for small magnetic fields. 


hence computing response functions involving the 
derivatives with respect to the magnetic field can 
be a challenging task (see for example Il20l in the 
context of magneto-electric response of topologi¬ 
cal insulators). With operator algebras, one has a 
well defined calculus 1211 , that is, a set of explicit 
rules which reduce these calculations to a simple 
exercise in differential geometry (see for example 
l22l for a computation of the magneto-electric re¬ 
sponse function done in this way, and to be used 
here). For the systems defined in |[2, these tools 
can be applied to investigate the behavior of var¬ 
ious quantities with respect to the parameter b |2 
(called 6 here). In particular, one can show that 
the invariants are continuous of b (hence constant) 
if a gap happens to remain open, a fact which was 
assumed in Q but questioned in ||7|, for example. 

A VIRTUAL CHERN INSULATOR IN 1+1-DIMENSIONS 
The System Defined 

The first example of a virtual topological insulator con¬ 
sists of a 1-dimensional periodic chain subjected to an 
incommensurate onsite potential: 

H^^Ti+Tl + 2Wcos{eX + (p), ( 1 ) 

where T„ represents the translation operator by n unit 
cells Tn\x} = |x + n), and X the position operator X|x) = 
x\x). The Hamiltonian has a localization-delocalization 
transition at the critical value Wc = 1: For W < 1, the en¬ 
tire spectrum is delocalized and for W > 1 the spectrum 
is localized [|23|. Throughout, we will set the fundamen¬ 
tal constants to one. 

This is the topological model considered in the Ref. 12 
and further analyzed in a number of subsequent works 
E HISHISI- While most of the related experimental 
works have been focused on photonics, here we adopt 
a different view and think of this model as that of a 
molecular chain absorbed on a substrate with incom¬ 
mensurate lattice. The angle 6 quantifies the substrate's 
lattice constant relative to that of the chain and the angle 
(p quantifies how the molecular structure is aligned with 
respect to the substrate's lattice. In a real experiment, 
one deals with an entire ensemble of absorbed molecu¬ 
lar chains, for which (p is expected to take random values 
in the interval [0, 27 t]. 

The bulk energy spectrum of as function of 6 is 
shown in the left panels of Figs. 1 and 2 for subcritical 
and overcritical values W = 0.5 and W = 1.5, respec¬ 
tively. The spectra are identical with the corresponding 
Hofstadter spectra (see discussion below). The energy 
spectrum of an ensemble of long but finite chains with 
open boundary conditions and random (j)-parameter is 
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FIG. 1. (Color online) Left panel: The energy spectrum of 
defined in Eq.j^as function of 6, computed with periodic 
boundary conditions on a finite 1 dimensional lattice of 10^ 
sites. Right panel: The energy spectrum of as function 
of 0, computed on a 400-sites lattice but with open boundary 
conditions. The plot shows a superposition of 100 spectra 
obtained with 100 random values of (p in the interval [0,2n]. 
For guidance, the bulk spectrum from the left panel is shown 
with darker color. The value of W is subcritical, W = 0.5. 


shown in the right panels of Figs. 1 and 2. They reveal 
that, regardless of fhe localized or delocalized characfer 
of fhe bulk sfafes, almosf all bulk energy gaps are com- 
plefely filled wifh edge specfrum. These are indeed fhe 
frades of a fopological insulafor. 


Connection with the Hofstadter Model 

Since Hp is nothing but the Harper or the almost- 
Mathieu Hamiltonian, its connection with models de¬ 
scribing two dimensional electrons in magnetic fields is 
well known. Here we provide an algebraic connecfion 
which will help us shape fhe principle for consfrucfing 
virtual topological insulators. Quite surprisingly, a key 
role is played by the angle cp. 

Note first that under translations: 

TnH^Tl = Ti + T| + 2 W cos(0X + (p + nd), (2) 

or: 

TnH,pTl - H^^+ngjmod271- (3) 

As a result, the translations generate a dynamical system 
(S, t) on the unit circle S: 

S 3 (p ^ T(p = {(p + 6) mod 2n. (4) 

Acting repeatedly with t, we obtain an action of Z on S: 

n times 

Z 3 n —> Tn(p =T o ... o T (p = {(p + n 6) mod 2n. (5) 

As we shall see, fhis d 5 mamical sysfem is key fo fhe self¬ 
averaging properfy menfioned in fhe infroducfion. 

Nexf, we define a dual d 5 mamical sysfem. For fhis, one 
considers fhe algebra C(S) of complex valued continuous 
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FIG. 2. (Color online) Same as Fig. 1 but for the overcritical 
value W = 1.5. 


functions over the unit circle. implements an action 
of fhe group Z on C(S) in fhe following way: 

C(S) 3/^ «„(/)=/ OT„. (6) 

Poinfwise, fhis is equivalenf fo: 

/((/)) ^/(T„(/)),V(^e [0,271]. (7) 

The dual dynamical sysfem is (C(S), a^. In the theory of 
operafor algebras, fhere is a sfandard consfrucfion which 
enlarges fhe algebra C(S) so fhaf fhe acfion defined in 
Eq.|^is implemenfed by unifary internal elements. It is 
called the crossed product of C(S) by Z and is denofed 
by C(S) Xa Z. If is obfained by augmenfing a unifary op¬ 
erafor u fo C(S) and enforcing fhe commufafion relation: 

w7w-” = «„(/), V / e C(S). (8) 

Then fhe elemenfs of C(S) ><« Z are formal series: 

a = Y^fnU\ fn^C{S). (9) 

If b - anofher element, then the commuta¬ 

tion relation in Eq. [pleads to: 

ub — ^ fm / (lH) 

«eZ meZ 

which should be seen as a convolution that is twisted by 
the map a. We should mention that the crossed product 
algebra C(S) Z is one of fhe mosf sfudied algebras in 
operafor algebra fheory ETI . 

The algebra C(S) Z accepfs a family of canonical 
represenfafions on fhe Hilberf space ^^(Z) of square- 
summable functions over Z: 

Ti^(a) = fn(T:x(p)\x}(x\T„, (pe[0,2n]. ( 11 ) 

fixXeZ 


If is a simple maffer of term counting for one fo convince 
himself fhaf indeed n(p{a)n^{b) - n^{ab). The family of 
Hamilfonians defined in Eq.j^are represenfations 

Hcj, = n^{h) 


( 12 ) 
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of the following element: 

h = u + u*+ 2Wcos{(p). (13) 


We now consider the Hofstadter Hamiltonian on the 


2-dimensional lattice Z^: 


H® = 7'(i,o) + + W^Tjo,!) + 1)) 

(14) 

where: 


7*(n,m)l^/ y) = + n,y + m) 

(15) 

are the dual magnetic translations and obey the follow¬ 
ing commutation relation: 

7'(i,o)‘7'(o,i) = c'®T(o,i)T(i,o). 

(16) 


In this form, the Hofstadter Hamiltonian is invariant 
with respect to the magnetic translations: 

+ n,y + m), (17) 


written here in the Landau gauge. Above, O represents 
the magnetic flux through the unit cell, expressed in units 
S/e. 

We show in the following that Ho is in fact a repre- 
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sentation on of the same element h from Eq. 

Indeed, since any continuous function over S can be de- 
cormosed as a discrete Fourier series, one can see from 
Eq.|^that the algebra C(S) ^ is generated by z = e“f’ and 
u. The commutation relation between these generators 


is: 


uz = e‘^zu, (18) 

which is exactly the commutation relation between the 
7~'s, if we enforce 0 = 0. If that is the case, we imme¬ 
diately obtain a representation of C(S) Xa Z on €{^), by 
sending z into 7~(i,o) and u into T(o,i). In other words: 

n'(a) = Y^fn(T(i,o))(r(o4. (19) 


It is evident that Ho = n'{h) with the same h as in Eq. 13 


Analysis of the Virtual Topological Insulator 

Given the cormection between H(f, and Ho, we can un¬ 
derstand at once why their energy spectra are identical. 
Indeed, recall that the abstract element h itself posses a 
spectrum, defined by those cj e C for which h - ■ 1 is 

not invertible. If 0 is irrational, then the representations 
71 and n' are faithful, and consequently the spectra of h, 
n^{h) and n'{h) coincide. We also know from Ref. 1281 
that the gap edges in the Hofstadter spectrum are con¬ 
tinuous with respect to 0, hence the spectra coincide 
also for the rational values of 0. One should be aware 


that this statement concerns strictly the location of the 
spectrum and that it says nothing about the localized or 
delocalized nature of the spectrum. 

The next thing we want to investigate is the nature 
of the virtual dimension and how to compute the bulk 
topological invariant. Using the generators z and u, we 
can write the generic elements as: 

fl = ^ CmnZ"'u"- (20) 

m,M€Z 

This representation reveals the 2-dimensional nature of 
the problem, as {m, n) live on the 2-dimensional lattice 
The coordinate n is related to the real-space coordi¬ 
nate X but the coordinate m is abstract and relates to the 
Fourier decomposition of the /„ functions. For example, 
the generator h of the Hamiltonian Hip is give by: 

h = u + u* + W(z + z*). (21) 

We now follow the work of Bellissard on aperiodic 
solids IITtII and introduce a non-commutative calculus on 
C(S) Z. It is defined by a formal (but mathematically 
rigorous) integration J: 

1(a) = Coo, (22) 

and two derivations: 

dia - i ^ mcmnz'"u" (23) 

m,n€.Z 

and 

02fl = t ^ nCmnZ"‘u". (24) 

m,n€Z 

In the representation given in Eq. the first derivation 
is just the ordinary derivation with respect to (p: 

dia = Y^{dipfn)u\ (25) 

neZ. 

and 

I(a)=ip^ fdcpMcp). (26) 

Js 

Together, (C(%) Z,0, define a non-commutative 

manifold dubbed m the non-commutative Brillouin 
torus. 

Now let y(t) = ^(1 + sgn(t)) be the usual step function 
on the real axis and consider the Fermi projection: 

p = x(eF - h), (27) 

that is, the spectral projection of h onto the spectrum 
below the Fermi level Cf. We define the bulk invariant 
as the 1-st non-commutative Chern number IfTSl : 

Chi(p) = 2nil{j>[dip,d2pi), 


(28) 
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which is a topological invariant 1181 in the sense that, as 
long as ep is located in a spectral gap of h, Chi(p) takes 
integer values and it remains constant under continuous 
deformations of h. 

As already elaborafed in fhe infroducfion, one of fhe 
main sfafemenfs of m was fhaf fhe bulk invariant can 
be computed at fixed (p. The implication of fhis sfafe- 
menf is fhaf fhe fopology is encoded in a single phys¬ 
ical sysfem rafher fhan in ifs whole family 
These sfafemenfs were reinforced in ||7l, in response fo 
fhe apparenfly confradicfing conclusions in ||6l- Now, 
examining Eq. 28 and fhe definition of fhe infegrafion 
in Eq. 26 fhere seems fo be a confradicfion here. Buf 


confrary we show fhaf fhe sfafemenf of HI is correcf and 
can be formulated in mathematically rigorous terms. To 
start, let us use Birkhoff's ergodic theorem ISJ to rewrite 
the integration as: 


ti|<N 


Then we use the representation TCcf, to write: 

fo(T:x(p) = {x\n^{a)\x} (30) 

and now one can see that the integration is nothing but 
the trace per length Tr^ of a physical represenfafion af 
fixed (p: 


^ ^ Zj = TrL(7i0(fl)). 


(31) 


|x|<N 


Eurfhermore, using fhe simple rules: Ti^(dia) = dipn,p{a) 
and n^{d 2 a) - i[X, ncj,{a)], we obfain fhe desired expres¬ 
sion for fhe bulk invarianf: 

Chi(p) = -271 Trp(P^[d^P^, [X,P^]]), (32) 

where P^ = xi^F - is fhe physical Eermi projection. 
The sfafemenf in HI is confirmed. A fine poinf which 
was made in iHllzl was fhaf fhe sfafemenf is frue only if 6 
is irrafional. We now can see explicifly fhe reason behind 
fhis sfafemenf: fhe map Zx is ergodic on S for 6 irrafional 
hence Birkhoff's ergodic fheorem applies (w.r.t. fhe usual 
measure), buf fhis is nof fhe case if 6 is rational. The 
former case can be visualized by imagining fhe sequence 
(p, T(p, T^cp,... being inked on a piece of paper, in which 
case one will discover fhaf, no maffer how sharp fhe 
writing insfrumenf, fhe unif circle will always be entirely 
covered by ink. This is obviously nof fhe case if 0 is 
rational, in which case fhe sequence of dofs will closed 
info ifself and fhere will only be a finife number of marks. 
Hence, fhe sysfem is a virtual topological insulator, as 
defined in fhe infroducfion, only if 6 is irrational. 


traduces random fluctuations in all the coefficients of fhe 
Hamilfonian: 

Ha,,4, = 2]^(1 + Aa;,)(|x)<x + 1| + |x)<x - 1|) (33) 

xeZ 

-f- 2W ^(1 + A'a;(.) cos(x 9 + (p) |x)(x|. 

xeZ, 

The disfribufions of ojx and oj^ should be defermined 
from fhe energetics of fhe lattice disforfions, buf here we 
will assume fhaf fhey are independenfly and randomly 
generafed from fhe inferval [-|, 5]. One should nofe 
fhaf fhe original Hamilfonian is recovered when A and 
A' are sef fo zero. 

As if follows from fhe works of Bellissard HHl, fhe 
disorder can be freafed algebraically foo. Eor fhis, one 
considers fhe disorder configuration space 

0 = [-l,lf x[-l,lf, (34) 

such fhaf CO = {co, co') e Q gives af once all fhe random 
amplifudes in fhe Hamilfonian: 

0.3 CO = \(Ox,co'x]x€-^- (35) 

The probabilify measure fo be used for disorder averag¬ 
ing is simply 


da) = ]^ (36) 

a: 

Nexf, one defines fhe ergodic d57namical sysfem (QxS, t), 
wifh: 


Tn{C0,(p) = {x„CO,Tn(p) 

= ({a;:c+„,a);+„};cez,((^ + n0)mod27i). (37) 

The relevance of fhis dynamical sysfem comes from fhe 
facf fhaf: 


T„Ha,,4,Tl = Hx^aj.intp- 


(38) 


Sysfems possessing such covariance properfy are called 
homogeneous. 

The dual d37namical sysfem (C(Q x S), is defined by: 

M/) =/°T„, V/eC(QxS). (39) 


Poinfwise, fhis is equivalenf fo: 

f{C0, cP) ^ f{T„CO, Tn(p), V (tU, (j)) G Q X S. (40) 


Robustness Against Disorder 

The disorder can occur in the chain's lattice or in the 
substrate. In a tight-binding approach, the disorder in- 


As before, one can form fhe crossed producf algebra 
C(Q X S) Xa Z by augmenting a unifary elemenf u fo fhe 
algebra of continuous functions over QxS, such fhaf: 

u>'fu-^ = an{f). 


(41) 
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The elements of this algebra are formal series of the form: Consider now the Fermi projection: 


a = J^fnu’\ /„eC(QxS), (42) 

n€Z 


v' = xiep -h')= Yj (50) 

m,ne'Z 


and, if 1; = another element, then the com¬ 

mutation relation in Eq.|^ leads to: 

Clb = ^ fm • ( 43 ) 

neZ meZ 

This looks formally the same as Eq. but note that f„ 
and a„ have a new meaning. 

The algebra C(Q x S) ><« Z accepts a family of canonical 
representations: 

no,,4,{a)= fn{TxCO,T^(p)\x}(x\T„. (44) 

n,j:eZ 


and it is easy to see that the disordered Hamiltonian in 
Eq. 33 is the representation 


bioj,(p ~ '^(o,cp{b ) 


(45) 


According to Ref. IITHl , as long as: 

Yj + n^) f dco \p'„n{a))\ < oo, (51) 

the first Chem number: 

Chi(p') = 2ni l(p'[d,p',d2p']) (52) 

remains quantized and constant under continuous de¬ 
formations of h'. This invariant can be evaluated nu¬ 
merically using the elementary methods developed in 
Refs. lf29lj3TI . Eor example, even in the regime where the 
Eermi level is imbedded in dense localized spectrum, the 
quantization of the Chern numbers can be obtained with 
machine precision when these algorithms are use. 

To understand the physical meaning of the mathemat¬ 
ical condition [M] first note that: 


of the following element: 


"^oj^cpip ) ~ (53) 


h' = {1 + Aa)o){u + «’) + 2W(1 -h A'co'q) cos{(p). (46) 

Here and throughout, by ojq we mean the function which 
assigns to m = {cOx]xez the real number coq- Note that 
when h' is multiplied by another element or, for example, 
is squared, the result of the multiplication (see Eq. |43) 
will contain = cOx, explicitly. 

Eollowing again the work on IQHE of Bellissard and 
his collaborators IITSl , we define a non-commutative cal¬ 
culus over the new algebra. Eirst, let us note that the 
elements of this algebra can be represented as in Eq. 
but this time the coefficients Cmn are not just numbers but 
functions over Q: 

a= Y^ (47) 

m,neZ 

The noncommutative Brillouin torus is defined as the 
non-commutative space 

(C(QxS) Xa Z,d,j) 

where the derivations remain the same as before but the 
integral is modified as: 

I{a) = [ dco cooico). (48) 

Jo 


i.e. the Eermi projection for one physical representation. 
Then: 


Pn{co,(p) = {0\P^^^\n}, 

and ondition|5T]can now be translated into: 
Y^n"^ f dcp f dco \{0\Pai,^\n)f 

v,tzl7 ^ 


< OO 


(54) 


(55) 


and 


Y^ f dcp f dco \d^(0\P^^^\n}\ 
«eZ y ^ 


< oo. 


(56) 


These two conditions must be simultaneously satisfied 
for the quantization and homotopy invariance of the 
first Chem number to hold. While Eq. 55 is automat¬ 


ically satisfied if the Eermi level resides in a region of 
Anderson-localized spectrum, this may not be the case 
for Eq. 56 However, the later condition is for sure satis¬ 
fied if the Eermi level resides in a clean spectral gap. 

Since the d 3 mamical system (QxS, t) is ergodic, we can 
follow exactly the same arguments as before to equiva¬ 
lently expressed the bulk invariant as: 

Chi(p') = -271 TrL{p^_p[d^P^^, [X,P^_^]]). (57) 


In the representation given in Eq. 42 

-^(«) = ii fd(p f dco fo{(o,(p). 
Js Jo 


As one can see, both integrals over <p and co disappear, 
which tells us that the invariant can be computed at 
fixed virtual coordinate cp and from one disorder config¬ 
uration. 
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Physical Prediction for Bulk 

The physical meaning of the invariant can be deter¬ 
mined by adapting the non-commutative Kubo-formula 
IITSl l32l l33l to the present context. This leads us to the 
prediction that by moving the chain relative to the sub¬ 
strate, hence varying the angle (p in time, one can set a 
charge current / along the chain. The time-average value 
of the current is: 

/ = Chi(p') ,p. (58) 

This is our physical prediction, concerning the bulk. This 
effect will be present in the absence of edges, such as in 
a ring geometry. 

Bulk-Boundary Correspondence 

We now consider the half-space system by restrict¬ 
ing the physical space from Z to N. Consequently, we 
replace the Hilbert space £^{Z) with ^^(N). As it fol¬ 
lows from the works of Kellendonk, Richter and Schulz- 
Baldes IITOlHill , this new physics can also be treated al¬ 
gebraically. To understand the required modifications, 
let us enquire about the faith of the translation operator, 
after we project onto the half of the space: 

Ti ^ ti = TiH, (59) 

where H is the projection from €^{Z) to f^(N). We have: 

ti|x) = |x + l), VxeN, (60) 

and 

t| \x) = |x - 1), V X > 0, and f\ |0) = 0. (61) 

As a consequence, the translation is no longer a unitary 
operator and instead is a partial isometry: 

t|ti=J, tit|=J-|0)(0|. (62) 

This suggests that the only modification we need to make 
to the algebra C(Q x S) ><« Z is to replace the unitary 
operator m by a partial isometry u: 

u*u - 1, uu* = 1 - e, (63) 

where e is a projection, = e, e* = e. This projection 
relates to |0)(0| once a physical representation is con¬ 
sidered. The elements of the new algebra, denoted by 
C(Q X S) Xa N (not a standard notation), are formal series: 

fl = fmn U^{UT, fmn ^ C(Q X S). (64) 

Now, any element a = Y^nez fn w” from the bulk alge¬ 
bra C(Q X S) Xa Z can be transformed into an element 
from the half-space algebra: 

a j{a) = 2], fn 

n€Z 


where by it is understood (m*)". Note that u is not in¬ 
vertible hence u~^ will be incorrect to use. Also note that 
j is only a linear map and not an algebra homomorphism 
(it does not respect the multiplication rule). Reciprocally, 
any element d - Ym.n^tifmn «"'(«*)" from the half-space 
algebra can be transformed into an element from the 
bulk-algebra: 

d ev(d) = Y fnn (66) 

m,n€N 

by sending u into u and u* into u~^. This time, the map ev 
respects both the addition and multiplication operations, 
hence it is an homomorphism of algebras. 

Using Eqs.|^ it follows immediately that any element 
from C(Q X S) x„ IN can be written uniquely as: 

d = i(a) + d, (67) 

with 

a = ev(fl) and d = d - i{a). (68) 

The component d is necessarily of the form: 

d= Y fn‘n HuT, fmn ^ C(Q X S), (69) 

m,M€N 

and the elements like d form an ideal inside the half¬ 
space algebra, which we call the edge-algebra. A formal 
integral can be defined over the edge algebra: 

I{d) = Y f f fnn{C0,(p). (70) 

„ Js Jo 

The elements of the half-space algebra accept a canon¬ 
ical representation as operators over ^^(N): 

d = i{a) + d-^ f[a,,(i,(d) 

^ nna,,4,{a)n + Y fmn{TmO),Tnt(p)fm\0){0\fl, 

m,nehi 

where na,,(p{a) is the bulk representation already dis¬ 
cussed above. We can see that the first term represents 
the restriction of the bulk operator na,,,p{,a) to the half¬ 
space via the open boundary condition, while the second 
term represents the component of the half-space opera¬ 
tor that is localized at the edge. For example, the family 
of half-space Hamiltonians is generated by: 

j{h') = (1 -t Aa)o){u + u) + 2W(1 + A'w^) cos{(p). (72) 

The open boundary condition can be changed to any 
other boundary condition by adding an element from 
the edge algebra: 

h=Y hmnU'’'£{uT. 

m,n€N 


(65) 


(73) 
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This term can account for the lattice distortions or re¬ 
laxations near the edge, chemical contamination, bond 
breaking and so on. The statements below are true for 
any such generic boundary condition. 

We now follow again Kellendonk, Richfer and Schulz- 
Baldes IITOlITT] , and define fhe fopological invarianf for 
fhe edge sfafes. We will consider fhe general case wifh 
disorder and an arbifrary boundary condifion af fhe 
edge. A bulk specfral gap is assumed. Hence, lef 
be a family of half-space Hamiltonians, generated by h' 
from fhe half-space algebra: 

h' = iih') + h', (74) 

wifh h' from fhe bulk algebra and given in Eq.|^ and h' a 
generic elemenf from fhe edge algebra. The consfrucfion 
of fhe edge invarianf sfarfs from fhe bulk Fermi projec- 
fion p' - xi^F ~ h'). The firsf fask is fo find an elemenf 
g from fhe half-space algebra such fhaf ev(g') = p'■ In 
mafhemafical ferms, we are searching for a liff of p' from 
fhe bulk fo fhe half-space algebra. In our contexf, fhis 
means: 


g'=j{p')+p'. (75) 

The solution is quite simple. Consider a smoofh ver¬ 
sion ic of fhe sfep function y, such fhaf y(e) = 1/0 
above/below a small interval around fhe origin (hence 
fhe enfire smoofh variafion happens inside fhe small in- 
ferval). Then g' can be faken as: 

g'=x(e,-h') (76) 

because, since fhe difference befween x x occurs 
inside fhe bulk energy gap, 

p' = x(eF - h') = m - h'). (77) 

and one can wrife: 

P' = m-h')-i(m-h')), (78) 

which finally can be shown fo be localized near fhe edge, 
hence belongs fo fhe edge algebra. This will nof be fhe 
case if x were replaced by x- In ofher words, fhe smoofh- 
ness of X is essential in this construction. 

The next step is to consider the unitary element: 

u' = exp ( - 2nig'^ (79) 

in the half-space algebra. Then fhe edge invarianf is 
simply fhe winding number of u' in one dimension: 

vi(m') = -l(u'd^{uy). (80) 

The firsf fundamental result of Ref. m is fhe equalify 
befween fhe bulk and edge invarianfs: 


(81) 


Lef us now consider fhe physical represenfafion: 

^CO,(p ~ )/ (^^) 

which can be computed direcfly from: 

Ua,,^ = exp ( - 2nix{eF - (83) 

by diagonalizing or by any ofher mefhods of func¬ 
tional calculus. Then the invariant takes a more familiar 
form: 

h(m) = -2S fd(p f (84) 

This is a frue edge invarianf because can be replace 
by - f in 84 Looking back af fhe consfrucfion, if 


immediately becomes apparent that Ua),cj) -1 can be built 
entirely from fhe specfrum and fhe sfafes inside fhe bulk 
energy gap, i.e. from fhe boundary sfafes. One can see at 
once that, if fhe bulk invariant is not zero, then the bulk 
spectral gap must be filled entirely wifh edge specfrum 
when (p is scarmed over fhe interval [0, 27 t]. Indeed, if 
fhere was a gap leff, fhen we can place fhe Fermi level in¬ 
side fhis gap in which case Ua,,^ -I - 0 and consequenfly 
vi(m') = 0. Buf fhis will confradicf fhe equalify befween 
fhe bulk and boundary invarianfs, hence no specfral gap 
can occur in fhe edge specfrum. Nofe fhaf fhis sfafemenf 
is now esfablished for any boundary condifion. 

We now show fhaf fhe boundary invarianf, hence also 
fhe bulk invarianf, can be refrieved from fhe boundary 
physics af a fixed (p and disorder configuration. An in¬ 
teresting remark in 13 was that if fhe boundary of one 
physical realization is efched afom by afom, then bound¬ 
ary states will appear and disappear in the process, and 
if we mark fhe energies of all fhese sfafes, fhen fhese 
marks will evenfually fill fhe insulafing gap enfirely. 
This will reveal indeed fhe fopological characfer buf if 
is nof clear how one will recover fhe bulk invarianf by 
jusf observing fhe edge eigenvalues (fhink of some large 
value, Chi = 10, for example). Now, when fhe afoms 
are efched, fhe boundary is moved buf we can shiff fhe 
whole sysfem (chain + subsfrafe) and place fhe bound¬ 
ary back in the original position. This set of moves will 
effectively shiff fhe disorder and advance fhe phase cp by 
d. Now consider fhe quanfify Tr{( Uco,(i, which 
can be computed enfirely from fhe boundary sfafes. If 
we average fhe proposed quantify as fhe boundary is 
efched afom by afom, fhen we can see again Birkhoff's 
ergodic fheorem in action: 


(85) 


x<N 

= -27 f ( 86 ) 

and fhe boundary invarianf emerges. The sfafemenf is 
proved. 


Chi(p') = vi(m'). 
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Physical Interpretation and Predictions 


The system defined 


Returning to physics, let us state the second funda¬ 
mental result of Refs. Goiini, which is fhe following 
equivalenf formula for fhe edge invarianf: 

(87) 


We consider fhe 4-dimensional rofafional algebra gen- 
erafed by fhe unifary elemenfs uq,. .. M 3 satisfying fhe 
commufafion relafions: 

UiUj = e'^‘'UjUi, dij = -dji e [0,27 t]. (91) 

This algebra is generafed by fhe monomials 


wifh p(e) - dgxie), f de p{e) = 1. Recall fhaf (p is acfually 
a disfance, hence 




( 88 ) 


is fhe force operafor. If p(Ha,^p) is inferprefed as a specfral 
weighf, fhen fhe above formula is nofhing else buf fhe 
quantized averaged boundary force acfing on fhe edge 
of fhe chain, discovered in Refs. lEllIl. 

One question fhaf arises is how is one going fo enforce 
fhe weighf disfribufion p on elections in practice? The 
answer is similar fo fhaf for IQHE: creafe an imbalance 
befween fhe occupation of fhe edge sfafes by applying 
volfage. Indeed, nofe fhaf if fhe chain has two edges, 
then at equilibrium there will be boundary forces on 
bofh edges which cancel fhemselves exacfly. Buf if one 
applies a pofenfial bias AV befween fhe fwo edges of 
a finife chain, fhen fhere will be an imbalance in fhe 
occupation of fhe edge sfafes and fhis creafes an effective 
weighf disfribufion on one of fhe edges: 


p(e) = 


lifee [ep- iAy,ef-h lAV] 
0 ofherwise. 


(89) 


We now can sfafe our physical prediction, namely, fhaf 
fhere will be a nef force on fhe chain, which, on average, 
is defermined by fhe bulk invarianf and by AV: 


(Fnet) 0 =Chi(p')Ay. (90) 

In fhe view of fhe above discussion, {Fnet}(p is also equal fo 
fhe average force exerfed sponfaneously af one edge, as 
fhe edge is efched atom by atom. The pofenfial difference 
AV can, in principle, be induced by irradiafing fhe chain 
wifh elecfromagnefic waves. 


A VIRTUAL CHERN INSULATOR IN 3+1-DIMENSIONS 


We now follow fhe program ouflined in fhe tirsf Chap¬ 
ter fo generate a virfual Chem insulafor in fhree physical 
dimensions and one virfual dimension. We primarily 
use fhis example fo show how fhe sfrafegy works in 
a new seffing. Below we describe fhe sfeps fhaf fake 
us from fhe absfracf seffing all fhe way fo fhe concrefe 
physical predictions. 


= ul\..uf, x = {xo,...,X 3 )eZ.\ (92) 

hence fhe elemenfs are formal series: 

a=YLfx M*/ (93) 

xeZ'i 

wifh fx jusf c-numbers. The crossed producf considered 
in fhe previous chapfer is isomorphic wifh fhe algebra 
generafed by uq and ui if we fake 6 = 601 ■ In facf, fhe 
whole rofafional algebra can be generafed as an iferafed 
crossed producf 1341 . 

The non-commufafive calculus is defined by fhe infe- 
grafion: 


J(«)=/o, (94) 


and by fhe derivafions: 

dju - i^^XjfxU^. (95) 

X€Z* 

Given a projecfion p in fhis algebra, fhe second non- 
commufafive Chem number is defined as IH: 

3 

Ch 2 (p) = A 2 2^(-l)‘'j(p n ^jp)' (96) 

crcSq j=0 

where S 4 is fhe group of permufafions and (- 1 )" is fhe 
sign of fhe perfurbafion a. A 2 is a proper normalizafion 
consfanf. If is known 1191 fhaf Ch 2 (p) fakes infeger values 
which are invarianf under continuous deformations of 
fhe projector p. 

The Xo'group of fhis algebra is known in any dimen¬ 
sion d (= 4 in our case). Kg = 'Z?‘‘ \ and if all difs are irra¬ 
tional, fhen fhe 8 generators of fhe Kg groups are known 
explicifly l34l . In particular. Kg is known fo confain a 
mulfifude of projecfions wifh non-zero second Chern 
number. Hence, fhe self-adjoinf elemenf: 

3 

h = J^{Uj + up + W(Mo + M^) (97) 

/=1 

will display a large number of specfral gaps and if we 
place fhe Fermi level in one of fhese gaps, fhen fhe Fermi 
projecfion p - xi^F ~ h) will have a non-zero second 
Chem number for many of fhese gaps. These affirma- 
fions can be easily tested numerically. 
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We now define a representation of the algebra on the 
Hilbert space which will lead us to the desired 

physical models. It is given by: 


Mo 


^ g/( 0 oX+,f.) 

n^{uj) =Tdi, i = 1,2,3, 


n^{uo) = e 


(98) 

(99) 


where 0o = (0oi, 0 O 2 , ^os), X = (Xi,X 2 ,X 3 ) is the position 
operator on and: 




( 100 ) 


with dj the generators of Let us verify the commu¬ 
tation relations between uq and m; for / > 1: 


e‘(®oX+9)7-^jy) ^ 


— g‘(OoX+(p)^-iY,k<d 
- eoi(yj+6,,}+ci))^-iZt<i dikVt^y + 


( 101 ) 


and 


Hence: 


- e'(Li=i So0j+ip)^-iLk<i8,kyk^y + 


Ti^{uo)n^{ui) = e'^°'n^(ui)nt^(uQ), 


( 102 ) 


(103) 


as required. The other commutations can be verified in 
the same way. The representation of a generic element 
Eq. 93 is given by: 


n^{a) = Y, 


(104) 


AreZ'i 


We now can write down our example of a Chem in¬ 
sulator in three physical dimensions and one virtual di¬ 


(105) 


mension. It is generated by h of Eq. 97 
Hcj, = n^(h). 

Explicitly: 

3 

H*=+ r].) + 2Wcos(0oX+ (/)). (106) 

/=i 

As one can immediately see, we are talking about a crys¬ 
tal in a magnetic field B which additionally experiences 
an incommensurate potential in all 3-directions. The 
magnetic flux through the unit cell's facet in the (jk) 
planne is djk in units H/e. If we consider the magnetic 
translations on written in the Landau gauge: 


T'^\y)=e'^^>i^>-^y + dj), 


(107) 


hence it is a homogeneous system. 

The bulk invariant for this system is given by the sec¬ 
ond Chern number in Eq. [^applied to the Eermi projec¬ 
tion p = x{£f ~ h). Our task now is to show that the bulk 
invariant is computable inside the physical dimensions. 
Eor this, let us write the invariant explicitly, using the 
physical representations n^{p) = xi^F - H^) = P,f,. Erom 


Eq. 104 


m=fo = 


d(b 

Jo ^ (109) 


hence: 


ch2(p)= a2 r S n 

aeSk '^0 j =0 


( 110 ) 


where doP^ = d^P^ and djP^ - i[Xj,P^] for j - 1,2,3. 
If doi angles are irrational, then the d 5 mamical system 
(p ^ (p + day is ergodic w.r.t. all three lattice translations 
and we can use Birkhoff's ergodic theorem 


-1 

dpfip) = y £/('/' + ®0l/) 

0 V€V 

to express the invariant as 

3 

ch2(p) = a2 Yi-'^r^P'P n 


( 111 ) 


( 112 ) 


CT€S4 


;=0 


where Try represents the trace per volume. This shows 
that the invariant can be indeed computed at fixed vir¬ 
tual coordinate, hence our model is a virtual topological 
insulator. 


Physical Interpretation and Predictions 

The bulk topological invariant is connected to the 
isotropic part of the magneto-electric response function: 


a - 


1 y ^ 


(113) 


where P is the macroscopic electric polarization vector. 
Indeed, according to Ref. Il22l which treated the generic 
case with arbitrary magnetic fields and aperiodic poten¬ 
tials, the change in a when p is varied by a whole cycle 
of 2n is equal to the second Chem number of the Eermi 
projection: 


An = Ch 2 (p). 


(114) 


then T', and Tdk commute and one can see that de- write this as: 


fined in Eq. 106 posses the covariant property: 


T;H^(T;r=H^^eoyk 


(108) 


pin 

\ f ‘ 

Jo 




= Ch2(p). 


(115) 
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Employing again Birkhoff's ergodic theorem and observ¬ 
ing that the integrand is a macroscopic quantity hence 
invariant to lattice translations, we arrived to: 


with fx now complex valued functions over Q. Each 
element accepts a representation as an operator on the 
Hilbert space 


3 o 2 p, 

d(pdBj 


2nCh2{p). 


(116) 


w»)=LE fx{TyC0) \y){y\ (125) 


This is our physical prediction. It can be tested experi¬ 
mentally in the following way Consider fhe small fime- 
modulafions: 


B(f) ^B + AB cos(2n/if) (117) 

(p{t) - (p + Ap cos{2nf2t). (118) 

Since fhe variafion of P wifh time gives fhe elecfric cur- 
renf ||35] : 


where n is fhe represenfafion already defined above. 
One can verify explicifly fhaf indeed noi,(p{ab) = 
Tioi,ctiici)na,,cf)ib), as required for a represenfafion. Eor ex¬ 
ample, fhe disordered Hamilfonian in Eq. 121| is gener- 
afed by 


h' = {1 + Amo) + wp -h (1 -h A'a)'o)W{uo + w^. (126) 

/=i 


7 = dtP, 


(119) 


The infegrafion becomes: 


one can see fhaf fhe two modulafions in Eq. 117 will 
induce AC charge currenfs af frequencies /i, / 2 , 2/i, 2 / 2 , 
/i + /2 and so on, buf af fhe frequencies f = fi ± f 2 the 
amplitude of fhe currenf is: 


= nHfi ± fi) Ch 2 (p)A<|,AB. ( 120 ) 


Robustness Against Disorder 

Let us now consider the effect of disorder, which in¬ 
duces random flucfuafions in fhe coupling consfanfs of 
fhe Hamilfonian: 

3 

h».=EE (1 AcOy)\y){y\{Td, + 7'J'.) 

j=l yal? 

+ 2wY,{ 1 + A'co'y) cos(0oj/ + p)\y){yl (121) 


I{a)= f dcofoico), (127) 

Ja 

and fhe derivafions remain unchanged. The second 
Chem number of a projecfion is formally defined as in 
Eq.|^ According fo Ref. 11911 , fhe second Chern number 
of fhe Eermi projecfion continues fo fake infeger values 
fhaf are invarianf under continuous deformations of fhe 
Hamilfonians, provided fhe Eermi resides in a specfral 
gap. As in fhe previous Chapfer, one can formulafe 
sfronger condifions fhaf cover fhe regime of sfrong dis¬ 
order. If can also be compufed af fixed virfual coordinafe 
and from a single disorder configuration. 

The non-commufafive second-Chern number can be 
evaluafed numerically using fhe mefhods developed in 
Refs. I 3 TI . Explicif numerical calculations can be found 
in Leung's PhD fhesis l3^ . 


The space Q of fhe disorder configurafion and fhe acfion 
Ty of on Q can be infroduced as before. The laffer is 
given by fhe shiff of fhe disorder configurafion by an y. 
The disordered Hamilfonian has fhe covarianf properfy 
wifh respecf fo fhe magnetic franslafions: 

( 122 ) 

The disorder is infegrafed in fhe algebraic approach 
by considering fhe algebra generafed by C(Q) and ufs, 
wifh fhe additional commufafion relations: 

f Uj = Uj if o Tdi), j = 1,2,3, (123) 

and /mq - wq// for all / e C(Q). The new algebra is 
generafed by fhe monomials f{co) and fhe elemenfs of 
fhe algebra are formal series: 


The Half-Space Algebra 

Lef us now consfrain fhe model fo half of fhe space, 
X 3 > 0. Eollowing fhe argumenfs from fhe previous 
chapfer, we define fhe half-space algebra as fhe alge¬ 
bra generafed by C(Q) and uq,. .. ,U 3 satisfying fhe same 
commufafion relations as before, wifh fhe excepfion fhaf 
M 3 is modified info a parfial isomefry: 

M 3 M 3 = 1 , U 3 UI = 1 - 63 . (128) 

The half-space algebra is generafed by fhe monomials: 

f{a))ul‘>upul^u’^{uiy = (129) 

wifh X = (xo,xi,X 2 ) e Z^. Hence fhe elemenfs are formal 
series: 


fl = 2^ fmn{C0,x)u^U^{u*3f. 

m,H€N xeZ? 




(124) 


(130) 
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Any element a - T^xez* from the bulk algebra 

can be transformed into an element from the half-space 
algebra: 

a i(a) = (131) 

{x,X3)€X* 

with u~^ interpreted as and reciprocally, any element 
a = fmn{M,x)u^u'^{ulY (132) 

m,n€N 

from the half-space algebra can be transformed into an 
element from the bulk-algebra: 

d —> ev(fl) = ^ ^ (133) 

m,H€N 

As before, any element of the half-space algebra can be 
written uniquely as: 

d = i{a) + fl, (134) 

with 


The Boundary Invariant 

We need to define a non-commutative differential cal¬ 
culus over the edge algebra in order to define the bound¬ 
ary invariant. This is given by the integration: 

m = Y f do} fn„{co,0), (142) 

«eN 

which is well defined only if is localized near the 
boundary, and by the derivations: 

djd = iY 2 ], ( 143 ) 

for i - 0 , 1 , 2 . 

We now restrict the bulk Hamiltonian to half of 
the space by imposing a generic boundary condition 
through an element h' from the boundary algebra: 

h'=i{h') + h'. (144) 


a - ev(d) and d-d - i(a). 


(135) As before, we define the element: 


The component d is necessarily of the form: 

d= Y 2],(136) 

m,«€N 

and the elements like d form an ideal inside the half¬ 
space algebra, which we call the boundary-algebra. 

The elements of the half-space algebra accepts a canon¬ 
ical representation as operators over x N), given 
by: 

fto,,cp{f) = n{Y f{TyCo)\y){y\)n (137) 


g' = Xi£f - h') (145) 

such that ev(^') = p', the Fermi projector of h', and then 
the unitary element 

u' - exp(-27if^'). (146) 

The boundary invariant is then defined as the 3- 
dimensional winding number of this unitary element: 

2 

V 3 (m') = A 3 2^(-l)‘'t( Yliu - l)dj(u' - 1)), (147) 

creSs j=0 


and: 

(138) 

na),<i,{uj) = TdjTl, j = 1,2,3, (139) 

with n the projection from to x N). The fam¬ 

ily of half-space Hamiltonians obtained from the bulk 
Ha),tpf via open boundary conditions, HH^y ^H are gener¬ 
ated by: 

3 

iih’) = (l+Aa)o)2]^(M; + i2*) + (l + A'a);)W(uo + Mo). (140) 

/=i 

The open boundary condition can be changed to any 
other boundary condition by adding an element from 
the edge algebra: 

h' = Y hmn{C0,x)u^u’^e{ulT. (141) 

m,M€N 

This term can account for the lattice distortions or re¬ 
laxations near the edge, chemical contamination, bond 
breaking and so on. 


with A 3 a proper normalization constant. According to 
Ref. 1371 , the non-commutative odd Chem number of 
unitary elements takes integer values and is invariant 
to continuous deformations. And according to Ref. ITTII 
(see Theorem AlO): 


V3(m') = Ch 2 (p')- (148) 


As in 1 + 1 dimension, the boundary invariant can be 
computed at fixed virtual coordinate by cleaving the sur¬ 
face layer by layer and averaging the quantity inside the 


formal integral in Eq. 147 


To determine the physical meaning of this invariant 
one needs a formula like Eq. 87 which is not available 
yet. Nevertheless, we point out that it is this boundary 
invariant that assures the existence of boundary states of 
the half-space Hamiltonian: 




(149) 


when (p is scarmed over the interval [0, 27i]. Indeed, 
the spectrum of h is the reunion of all spectra of 
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and if this spectrum is gapped, then we can deform x 
such that its variation occurs entirely inside this gap. 
Consequently, up - 1 = 0 and the boundary invariant is 
zero. Then equality |148| ensures that, whenever the bulk 
invariant is no zero, the bulk gap is entirely filled with 
boundary spectrum. 

CONCLUSIONS 

Our conclusion is that Refs. |1T1|2| have indeed intro¬ 
duced a new class of strong topological systems. We are 
particularly excited about it because we are no longer 
bound to the 3 physical dimensions, which will certainly 
lead to the discovery of many new strong topological in¬ 
sulators. It will be interesting to explore the other sym¬ 
metry classes and see exactly how the virtual topological 
insulators are constructed there. 
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